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ABSTRACT 

The X-ray intensity of knot HST-1, 0.85" from the nucleus of the radio galaxy M87, has increased by more 
than a factor of 50 during the last 5 years. The optical increase is similar and our more limited radio data 
indicate a commensurate activity. We give the primary results of our Chandra X-ray Observatory monitoring 
program and consider some of the implications of this extreme variability in a relativistic jet. We find that the 
data support a 'modest beaming synchrotron' model as indicated in our earlier papers. Based on this model, the 
decay of the X-ray lightcurve appears to be dominated by light travel time across the emitting region of HST-1, 
rather than synchrotron loss timescales. 

Subject headings: galaxies: active — galaxies: individual(M87) — galaxies: jets — X-rays: general 



1. INTRODUCTION 

This paper is the third from a project to monitor the M87 
jet with the Chandra X-ray Observatory with complimentary 
data from the Hubble Space Telescope (HST), and more re- 
cently with the Very Large Array (VLA) and the Very Long 
Baseline Array (VLBA)^. The main goal of the project was 
based on the premise that synchrotron X-ray emission would 
come from electrons with halflives of order a year or less 
whereas inverse Compton X-ray emission would come from 
very low energy electro ns with halflives exceeding 10"^ years 
jHan-is & Krawczvns ki"2002'). In 'HaiT is et all (12003) (here- 
inafter 'Paper F) we reported the results from the first year, 
2002Jan-Jul. The main findings were significant intensity 
variability of the core (expected), and of HST-1, a knot in 
the jet 0.85" (65pc projected) from the core. We argued that 
the intensity changes were larger for the harder X-ray bands 
and suggested a synchrotron loss model with modest beaming 
(Doppler factor, S, of order 3 or 4, in agreement with earlier 
proper mot ion studies: piretta. Sparks. & Macchetta 11999) ). 
In Perlma n et alJ (Paper II), it was suggested that the 
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timescale for optical decay of the lightcurve for HST-1 was 
similar to the X-ray decay time, and thus not the factor of 10 
longer which was predicted by our original synchrotron loss 
model. Since the measured decay time was based on only a 
single minor change of optical intensity for each color, we 
consider this question still open. 

In this paper we concentrate on the recent Chandra data 
from 2003-2005 but will also include optical and radio data 
for HST-1 as required, leaving the details to other papers in 
preparation. The 'giant flare' of HST-1 (the X-ray intensity 
increased by a factor of 50 in 4 years) has forced us to abandon 
our original methods of data reduction which were devised the 
first year of our monitoring when ACIS pileup was not a prob- 
lem. ['ACIS' is the acronym for the Advanced C CD I maging 
Spectrometer and 'pileup' is discussed in section]^]] After 
a short section on the observations (§ 2) we deal with the data 
reduction methods and give the basic results for each type of 
measurement (§ 3). We review the synchrotron parameters in 
§ 4, and discuss implications of the lightcurve timescales in 
§5. 

We take the distance to M87 to be 16 Mpc ( lTonrvll991l) so 
that one arcsec = 77pc. We use the conventional definition of 
spectral index, a: flux density, S,y oc v'". 

2. OBSERVATIONS 

The first Chandra observations of M87 that specifically 
targeted the jet were those done by Wilson in 2000 Jul 



2 



Harris et al. 




Fig. 1. — The inner section of tlie M87 X-ray Jet: a staclced image from epochs A to V with an effective exposure time of order llOks. The data come from 
evt2 files so the true (average) intensity of HST- 1 is underestimated because of pileup. The coordinates labels are Right Ascension and Declination in J2000. The 
intensity units are electron volts per second per 0.0492 arcsec pixel, and the peak value is 5.5. The faintest red corresponds to 0.5 and the black below the green 
is close to 0.01. The yellow rectangle around HST-1 is that used for the lightcurve photometry and the small circle is that used in Papl (sections 3.2, 3.3.). 



JWilson & Yang|l200l iPerlman & WsodII2001 . They ex- 
perimented with various ACTS frame times in order to mini- 
mize the effects of pileup (two or more photons arriving at the 
same pixel during the same frame time, thus masquerading as 
a single event with the sum of the energies; see, e.g. Davis 
ll2001l) '). They found that a frame time of 0.4s gave essen- 
tially the same countrate as 0. 1 s for the nucleus of M87 which, 
at that time, was the brightest feature. Thus we adopted the 
same ACIS setup: chip S3 (back illuminated for low energy 
response), l/8th subarray, and 0.4s frame time. A list of our 
observations is given in Table[2 Basically, we have been get- 
ting a 5 ks observation approximately every 6 weeks for four 
'seasons': mid November to early August (when M87 is far 
enough from the Sun to be observable with Chandra and the 
HST). Our approved proposal for Chandra A07 will continue 
this strategy (2005Nov to 2006Aug). During the 2005 season 
we have obtained an extra 5 observations designed to give us 
weekly coverage for 2 months to check for shorter timescale 
changes. 

3. DATA REDUCTION AND BASIC RESULTS 

While we have maintained a fixed observing strategy, the 
data reduction has been done over the course of several years 
and thus the software reduction package, 'CIAO'^ (Chandra 
Interactive Analysis of Observations), has evolved through 
several releases, as has the calibration database. However, 
all data except the 'keV/s' values have been corrected for the 
contamination on the ACIS detector, and we are not aware 
of any software or calibration changes that would adversely 
affect our results. 

The essential preparatory procedures are removal of the 
pixel randomization and registration to align the X-ray core 
with the radio core. After setting up the bad pixel file and 
checking the lightcurve of a large background region for flares 
(since we use small extraction regions and no large flares were 

^ |http : //asc ■ harvard ■ edu/ciao/| 



found, we did not time filter any observation), we removed the 
pixel randomization imposed by the standard pipeline by fol- 
lowing the relevant CIAO thread. Registration of the X-ray 
image was accomplished by measuring the apparent positions 
of the core and knot A on both radio and X-ray images (the 
latter binned to 1/1 0th native ACIS pixel and smoothed with 
a Gaussian of FWHM=0.25"), and then changing the values 
of the FITS header keywords in the Chandra event file. In 
general, this procedure introduces an uncertainty of approxi- 
mately 0.05" or less for the early part of the program when 
the intensity of the core was comparable to that of HST-1. 
However, when HST-1 became more than 10 times brighter 
than the core, our ability to determine the centroid of the core 
emission is compromised, and we then gave more weight to 
the position of knot A which is not so well determined be- 
cause knot A is resolved. 

3.1. Intensities from fluxmaps 

The original data reduction procedure is outlined in Paper I 
and involves constructing monochromatic exposure maps for 
three energy bands. Dividing the data by the exposure map 
we obtain 3 flux maps. Each flux map was then multiplied by 
hvjtom to obtain units of ergs cm"^ s"'. v^q^ corresponds to 
the nominal energy for each band (chosen to approximate the 
band center weighted by the effective area) which was used to 
construct the relevant exposure map. This permits us to mea- 
sure fluxes with normal photometric procedures. All that is 
required is to sum the flux within a measuring aperture and 
multiply the result by the ratio of the mean energy of events 
within the aperture to the nominal energy of the band being 
measured. We have redone the fluxmaps now that CIAO au- 
tomatically corrects the exposure maps for ACIS contamina- 
tion buildup. Previously we had used a fixed region of M87's 
thermal hot gas emission to make this correction. 

When using small apertures (our standard region was a cir- 
cle with r=0.44" in order to separate the core from HST- 1 ; see 
fig. [3 a point spread function (PSF) correction of 1.55±0.05 
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TABLE 1 

Chandra dates and X-ray intensities of HST- 1 





Observational Parameters 


Intensities of HST-1 


Epoch Date 


Obsid 


Livetime 


K" 


f(0.2-6keV)'' 


S(2keV)'^ 


label 






(sec) 


(keV/s) 


CIO-'- ess) 




A"! 


2nnnTiinn 


1808 


12845 


0.247 


0.729 


29.2 




ZUUiJ all I yj 


3085 


4889 


0.734 


2.32 


92.8 


Q 


9nn9Fphl9 


3084 


4655 


0.588 


1.87 


74.8 




z.uwZ.ividJ. jw 


3086 


5089 


0.576 


1.84 


73.6 






3087 


4973 


0.798 


2.57 


102.8 


p 


90091,, 19/1 


3088 


4708 


1.096 


3.54 


141.6 


a 


9nn9i\rr\v 1 1 


3975 


5287 


799 


2 62 


104 8 


JJ 




3976 


4792 


0.653 


2.14 


85.6 


J 




3977 


5276 


0.645 


2.12 


84.8 


T 
J 


zXj\Jji\ld.l\Jy 


3978 


4852 


0.872 


2.88 


1 15.2 


]^ 


900"? Anrl A 


3979 


4492 


1.071 


3.54 


141.6 


L 


9nn^Mavl s 

Z.\JU-.'iVla.y 1 o 


3980 


4788 


1.022 


3.39 


135.6 


M 


9nn^Tiiin^ 

Z.U\JJJ \llUj 


3981 


4677 


0.859 


2.86 


1 14.4 


N 


9onsAiipns 

^ UV/ ^ U. t v/ o 


3982 


4841 


1.214 


4.06 


162.4 





2003Novll 


4917 


5028 


2.192 


7.36 


294.4 


p 


2003Dec29 


4918 


4677 


2.041 


6.88 


275.2 


o 


2004Febl2 


4919 


4703 


4.079 


13.81 


552.4 


Re 


2004Mar29 


4920 


5235 








s 


2004Mayl3 


4921 


5251 


5.358 


18.23 


729.2 


T 


2nn4Tiin23 


4922 


4543 


5.323 


18.16 


726.4 


u 


9nn4Aiion's 

Z. U W+rt. Ll g U J 


4923 


4633 


5.636 


19.28 


771.2 


V 


2004Nov26 


5737 


4237 


7.494 


25.76 


1030 


w 


2005Jan24 


5738 


4666 


8.316 


28.66 


1146 


X 


2005Febl4 


5739 


5154 


8.785 


30.30 


1212 


Ya 


2005Apr22 


5740 


4699 


12.42 


42.91 


1716 


Yb 


2005Apr28 


5744 


4699 


12.17 


42.08 


1683 


Yc 


2005May04 


5745 


4705 


11.80 


40.81 


1632 


Yd 


2005Mayl3 


5746 


5142 


11.60 


40.14 


1606 


Ye 


2005May22 


5747 


4701 


11.24 


38.92 


1557 


Yf 


2005May30 


5748 


4699 


10.66 


36.92 


1477 


Yg 


2005Jun03 


5741 


4698 


10.43 


36.14 


1446 


Yh 


2005Jun21 


5742 


4703 


10.27 


35.58 


1423 


Yi 


2005Aug06 


5743 


4672 


7.105 


24.64 


985.5 



Note. - All observations had l/8th subarray ACIS-S7 chip only, and 0.4s 
frame time. 

" The values of K are our primary measured quantity. They come from sum- 
ming the energies of all events (from the 'evtl' file) within a rectangle of 
length 5 pixels transverse to the jet and 2.5 pixels along the jet (fig.0^ The 
upstream edge is approximately midway between HST-1 and the core com- 
ponent. No background subtraction was employed: when the source is weak, 
the background (in counts) is of order 1 %, and for the highest intensity, the 
background is less than 0.3%. There is no correction for the buildup of ACIS 
contamination. 

f(0.2-6keV) is the observed X-ray flux (the flux density integrated from 0.2 
to 6 keV) with units of 10"'^ erg cm"^ s"'. An approximate coiTection has 
been made for ACIS contamination buildup but not for Galactic absorption. 

The flux density at 2 keV has been calculated from the observed flux assum- 
ing a spectral index, Q!=1.5. For «=1.75, the listed values should be reduced 
by a factor 0.75 and for a=1.25, they should be increased by a factor of 1.25. 



Hz- 



I njy is 10 

Archival data from lWilson & Yang|l20Ql . 
° This observation was taken in continuous clocking mode, so there is no 2D 
image available. 



is required to 'recover' the fraction of the PSF outside the 
small circle. Although it is known that the PSF size is a func- 
tion of energy, from our PSF simulation with ChaRT/Marx, 
we found values between 1.50 and 1.60 for our 3 energy 
bands. From the same simulations, we find that 6±4% of 
each PSF falls within the adjacent measuring aperture (this is 
for the special case of the small circles for the core and HST- 
1). In Paper I, these corrections had not yet been devised, so 



the flux scales in the figures should not be used to derive ab- 
solute values. We provide our best flux estimates later in this 
paper. 

3.2. Intensities from piled up data: keV/s 

Pileup for CCD's occurs when two or more events are re- 
ceived for the same pixel during a single frame time ( Davi^ 
1200 Ih . Thus we lose the inherent spectral energy distribu- 
tion one normally associates with ACIS data. However, to 
first order, the energy content of the arriving photons is pre- 
served and a measure of the total band is still available. To 
obtain intensities for lightcurves, we first followed the 'after- 
glow thread"* to recover the events with grade migration that 
caused them to be mistaken for afterglow events. Then we 
checked to see that all events of interest would be accepted by 
a suitable energy filter The incident spectrum of HST-1, both 
as measured when the knot was weak (essentially no pileup) 
and from our observation in continuous clocking mode, is es- 
sentially contained in the band 0.2 to 6 keV (probably 95% of 
the counts are below 4 keV). To recover the energy associated 
with piled events that now appear above 6 keV, we examined 
the event file at the higher energies and set the upper cutoff to 
a value that includes all events associated with HST-1 (up to 
17 keV in some cases). 

The CIAO tool 'dmstat' then gives the total energy (in eV) 
of all events. Dividing this by lOOOxlivetime (where iive- 
time' is the total time the detector is accumulating data) pro- 
vides an intensity measurement in ke V/s. Through early 2004, 
pileup was fairly mild, so we continued measuring the event2 
file. As the pileup increased, we found we were losing more 
and more events to grade migration so the data shown in fig.|2l 
are based on the event 1 file with no grade filtering and only 
the standard good time intervals applied. A rectangular ex- 
traction region (see fig-^* was used for this photometry. 

There are several limitations to this method. It is impossible 
to get spectral information like band fluxes; there is no a pri- 
ori method to correct for ACIS contamination buildup^; and 
there may have been bona fide events that were rejected by on- 
line software aboard the satellite and thus are unrecoverable. 
In spite of these deficiencies, this 'keV/s' method provides a 
reasonably direct estimate of the total observed intensity. 

3.3. Problems Associated with the Spectral Analyses 

For spectral analyses of HST-1 there are between several 
hundred and more than 3000 counts in a typical 5ks obser- 
vation. We have performed spectral fits for simple power 
laws (PL), power laws with a high energy cutoff (cutPL), 
and both mo dels in conjunction with the 'jdpileup model' 
(iDavisI 120011) . as implemented in Sherpa (.Freeman et alJ 
2001). There are a number of problems. 

The jdpileup model, even with the 'Monte-PowelF method 
which is advertised as being the best option to locate the true 
minimum in chi-squared, can produce a best fit which is ob- 
viously wrong: the fit gives an amplitude which is up to a 
factor of 10 too large, with an excessively large pileup frac- 
tion. This formal fit occurs because the fit finds a solution on 
the heavily piled up section of the pileup curve ( Davis 2001] 
fig. 6). This is particularly evident for featureless spectra such 
as power laws. In the fitting algorithm, there is a 'degeneracy' 

http : / /cxc ■ harvard ■ edu/ciao/threads/acisdetect afte rglow/1 
Note however that with the advent of pipeline version DS7.4, this procedure 
is no longer required. 

^ http : //cxc ■ harvard. edu/ciao/ threads /aciscontam/l 
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Fig. 2. — Multi-wavelength lightcurves for knot HST-1 in the M87 jet. The left panel shows the X-ray, and UV lightcurves on a linear scale to more clearly 
show the higher levels. On the light is the same plot on a log scale in order to better compare slopes (percentage changes). For this version we have added radio 
data at 15 GHz (the 2000 datum comes from an archival VLA observation). Both X-ray curves are from the evtl files in keV/s. For the UV data, the left panel 
has units of 0. 1 mjy whereas on the right, the units are mJy. 



for featureless spectra between the power law exponent and 
the pileup 'survival probability' parameter, jdp. alpha. 

The only way we found to mitigate this problem was 
to constrain several fit parameters: Nh, the total column 
density of neutral hydrogen per square cm towards the 
source, was for ced to be at lea st the Galactic value of 
2.5 xlO^" cm"2 JStarket al.lll992h : the survival probabiHty 
parameter jdp. alpha.min=0. 2; and the spectral index, a, was 
constrained to lie between and 3. It was also necessary to 
pick good values for the initial guess for jdp. alpha, a, and the 
amplitude of the power law. 

Another uncertainty is the mutual dependence between Nh 
and a which can be demonstrated by a simple example us- 
ing the spectral file for HST-1 from epoch D. The power law 
fit gives Nh=(3.4, H-1.9, -1.5)x lO^'W"- and a=1.46±0.14. 
Freezing Nh to values of 2.5, 3.5, and 4.5 x 10^*^, results in a 
values of 1.39±0.07, 1.47 ±0.07, and 1.54±0.08 respectively. 
Although it is conceivable that Nh could vary because HST- 
1 is at a projected distance of only 65pc from the nucleus of 
M87 and its diameter is less than a parsec, we assume it is con- 
stant, possibly augmented by an excess absorption above the 
Galactic value, i.e. intrinsic to M87. From spectral fits before 
pileup became a problem, we find fit values of Nh which av- 
erage 3.5 X 10^'^cm"^ (i.e. a bit larger than the Galactic value), 
although each individual solution is consistent, within the er- 
rors, with the Galactic value. 

A further effect complicates the solution. The pileup effect 
we call 'Eat Thy Neighbour' comes into play for two close 
sources. If a core photon arrives during the same frame time 
as a photon from HST-1, often they will be within each other's 
3x3 pixel region and thus be interpreted as a single event at the 
location of the pixel with the larger energy. 

Finally, the value of a depends on the choice between a 
simple power law and a power law with an exponential cut- 
off. Although in many cases the value of the cutoff energy is 
very large and the a's of PL and cutPL are sensibly the same, 
for other cases for which the fit value of the cutoff energy is 



below 12 keV, the corresponding values of a for the two mod- 
els diverge with a concomitant shift in the fit value of Nh (an 
example is given in Table |2}. 

Considering these difficulties, we believe it would be pre- 
mature to present a table of spectral parameters. For the 2002 
to 2003 data (prior to the large intensity increase), with Sherpa 
we find a «1.5±0.3, and Nh appears to be consistent with, or 
a bit larger than, the Galactic value. A joint fit of these data 
with XSPEC constraining Nh to a common value but allowing 
amplitude and spectral index to vary for each observation, re- 
sults in Nh=(2.6±0.5)x 10^" cm"^ and a ranging from a min- 
imum of 1.1 to a maximum of 1.5. It is our opinion that with 
spectral fitting, we cannot confidently claim that the spectrum 
of HST-1 is harder or softer when the intensity is high. 

Given the sizable uncertainty in spectral parameters of 
HST- 1 from piled up data, we proposed a 1 Oks observation for 
Director's Discretionary Time to obtain a continuous clocking 
(CC) mode observation to be made during those few days in 
2004March when the CCD readout direction is orthogonal to 
the jet. The idea was that we could obtain one set of reliable 
spectral parameters and compare them to those obtained from 
our normal 5 ks observation which was already scheduled for 
the last week of March. Our proposal was not granted time, 
but we were allowed to change the mode of our normal 5 ks 



TABLE 2 

Spectral Parameters for the CC mode data on HST-1 



Model 


Nh 


OL 






1020 cm-2 






power law (PL) 


7.3±0.9 


1.72±0.06 


0.42 


PL (BGsubtract) 


6.6±0.9 


1.73±0.06 


0.40 


PL + cutoff 


2.5±0.4 


1.18±0.11 


0.51 



Note. - For the PL -I- cutoff, the high energy cut value is 3.7itl.l keV and 
the value of Nh is the lowest allowed (set to the Galactic value). 
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observation to CC mode. 

In CC mode pileup is not a problem since the readout is con- 
tinuous and each pixel is effectively moved every 41 /^s. We 
followed the thread to correct the time of each event from the 
readout time to the arrival time, thus ensuring that the correct 
aspect solution would be applied. A profile along an East- 
West line (i.e. the jet is projected onto the RA axis from its 
normal position angle 20° away) is shown in fig.|3l We then 
isolated knot A, HST-1, and the core with thin rectangles on 
the ID 'image' of the event file. Notwithstanding the fact 
that there are no accurate calibration files for CC mode, we 
generated 'pha files' (the collected events with their energies) 
and corresponding effective area and redistribution matrices 
to permit us to run normal Sherpa spectral analysis. Back- 
ground areas on each end of the jet were used. The spectral 
results are given in Table |2l 



TABLE 3 
Doubling Times for HST-1 



Epoch Interval 


At 


(y-if 


dt" 








(yrs) 


E-F 


0.125 


0.373 


0.34 


I-J 


0.090 


0.352 


0.26 


M-N 


0.099 


0.413 


0.24 


P-Q 


0.123 


0.998 


0.12 


X-Ya 


0.183 


0.413 


0.44 



Note. - With a statistical error of a few percent on each intensity measure- 
ment, we estimate the uncertainties of the doubling times are dominated by 
systematic errors of order 15%. 
" y is the intensity ratio, hfli . 
dt is the doubling time in our frame. 




3994.050 4004.050 4014.050 4024.050 

4082.853 4082.555 4082.258 4081.960 

physical coordinates 

Fig. 3. — The profile of the jet along the RA axis of the observation of 
2004Mai'29 (Epoch R) taken in the continuous clocking mode. The resolu- 
tion is 0.1 native ACIS pixels, or 0.0492". Emission from the core, with an 
intensity similar to knot A (the feature towards the right), is lost on the left 
flank oi^ HST- 1 , the dominant feature. The location of the core should be 1 .6 
pixels to the left of HST- 1 . Note the relatively poor s/n of knot A, since the 
profile contains the background from the whole chip. 

The CC mode observation was meant to confirm the jdp 
spectral parameters, but insofar as a is concerned, the CC fit 
gives a significantly larger value than the preceeding and fol- 
lowing values. At this time, we do not know how to decide if 
the CC value is wrong because of unreliable calibration files, 
or if the jdp values are (systematically) too small. We note that 
whereas for the normal 2D image analysis, the background 
level has essentially no effect on the spectral fit since the mea- 
suring aperture is so small, this is not the case for the CC 
mode, for which the background of each pixel comes from an 
entire column in chip coordinates. While we cannot be certain 
that the spectrum of HST-1 did not change markedly during 
the months surrounding the CC observation, this seems un- 
likely. Although the CC experiment failed to provide higher 
accuracy for the spectral parameters, the range of values in 
Table|2]is still within the ranges for HST-1 prior to the major 
intensity increase. 

3.4. Fluxes and flux densities 

While most fluxes quoted from Chandra data are based on 
the usual spectral analyses (sec. 13. 3> . we favor fluxes from 
fluxmaps as being a quantity which is closer to an 'observ- 
able' and less susceptable to model assumptions. As de- 
scribed above, measuring on flux maps achieves this objective 



when pileup is not a problem. To generate approximate fluxes 
from the observed keV/s measurements, we convert keV to 
ergs and divide by the effective area. 

Since all indications are that most of the flux falls within 
the 0.2 to 6 keV band, we take a nominal effective area of 400 
cm^ as a first guess for the spectrum of HST-1 weighting the 
Chandra effective area. We then calculate a flux: 

F(0.2 - 6 keV) = 1.602 x lO-'^^^fri = 4.00 x 

JQ-I2 I gj.g ^^-2 g-I 

where K is the measured keV/s value and a is a correction 
factor for the (unknown) effective area which is a function of 
time because of the ACIS contamination buildup. 

We find the values of a by comparing F(a=l) with the 
flux measured on the flux maps for those observations made 
in 2002 and 2003 before pileup became a serious problem. 
We find a(2000.5828)=1.356, thereafter decHning to values 
around 1.2 by mid 2003. Making a rough extrapolation of a(t) 
permits us to make a first order correction for the contamina- 
tion problem. We estimate that resultant uncertainties are of 
order 15%, based on the scatter of a(t) from a smooth curve. 
Fluxes are listed in Table ^ In the same table we give flux 
densities at 2 keV calculated from the flux values assuming 
a=1.5. 

4. SYNCHROTRON PARAMETERS 

In this section we examine the physical parameters associ- 
ated with the synchrotron emission model in the same context 
as used in Paper I: we use the rise time to estimate the physical 
diameter and the decay time to constrain the beaming factor, 
assuming losses dominate. 

For a given interval. At, we measure a fractional increase 
of y=l2/Ii. Then the doubling time, dt= [(l/(y -l)]xAt, and 
the doubling time in the jet frame is dt' = Jxdt yr. Thus the 
diameter of the emitting region will be d < dt'/3.26 pc and 
the radius, r = d/(2x77) arcsec. From inspection of fig. |2j 
we select a number of intervals with large rate of increase 
and in Table |3l we give the calculated doubling times. For 
a selection of values for the beaming factor, in Table |4] we 
give the corresponding upper limit to the sizes of the emitting 
region for a characteristic doubling time of 0. 14 yr. 

We construct a spectrum for 2005 Jan 1 with data from the 
VLA; at 220nm from the Hubble Space Telescope; and our 
Chandra value for Jan 2005 from Table ^ This is shown in 
fig.|4] We then convert the frequencies and flux densities to 
the jet frame and solve for the usual synchrotron parameters 
JPacholczvkll97(]h : see TableEl 
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10= 10'» 10" 10'= 10'= 10" lO's 10i« 10" 10" 
Frequency (Hz) 

Fig. 4. — The spectrum of HST-1 around 1 Jan 2005. The radio flux 
densities are from the VLA; the UV from the HST; and the Chandra interpo- 
lated from Table m The two point spectral indices are «,-_uv=0.45 (22GHz 
to 220nm) and aox=1.03 (220nm to 2 keV). 



We have identified the most rapid decay rates, and these are 
given in Table|5] Estimates of the time for the intensity to drop 
by a factor of two range from 0.2 to 0.3 years. Matching 0.2 
yrs to the resuhs of Table|5j we find 6=6.3, a value somewhat 
higher than that found in Paper I and also larger than the mod- 
els reported in Table 3 of Biretta. Sparks. & Macchetto ( 199^ 
for the 6c proper motions in HST-1. However, the X-ray de- 
cay time is remarkably close to the rise time (which we used to 
estimate the source size), so it seems likely that both rise and 
fall reflect the time required to cross the source, and not the 
half-life. Although this is not proven, for the remainder of this 
paper we adopt the view that the X-ray decay time reflects the 
light travel time across HST-1. Thus solutions for smaller S 
(e.g. 3 to 5) are favored. Given the assumptions of the model, 
Tj'=diameter(l.y.) for 6 « 4.6 (see the next to last column in 
Table |5}- For smaller 6 the lightcurve decay should be dom- 
inated by light travel time; for larger S the weaker magnetic 
field leads to longer E" halflives. Large 6 (e.g. 20) values 
are disfavored because of their large r values and the much 
smaller angles to the line of sight which would be required. 

In Paper I we argued that the total emission of HST-1 did 
not impose any significant drain on the total estimated energy 
flow down the jet. This continues to be the case so long as 
(5 is > 3. The power needed to generate the required total 
energy in the time dt' (the jet frame doubling time) is listed 



TABLE 4 

Size of Emitting Region for HST-1 



S"= 


1 


3 


5 


9 


12 


20 


40 


dt"'(yrs) = 


0.14 


0.42 


0.70 


1.3 


1.7 


2.8 


5.6 


diameter(pc) = 


0.043 


0.13 


0.22 


0.39 


0.52 


0.86 


1.72 


radius(arcsec) = 


0.00028 


0.00084 


0.0014 


0.0025 


0.0034 


0.0056 


0.011 



Note. - For an observed doubling time of faO. 14 years, we calculate 
maximum values for the source size. 



"5 is the beaming factor. 

''dt' is the doubling time in the jet frame. 



TABLE 5 

Synchrotron Parameters for HST- 1 



s 


radius^ 


Deq 




log Eto," 




~uv 


r,'/d(l.y.)*! 






(cm) 


(gauss) 




(ergs) 


(erg s"' ) 


(yrs) 


(ly) 


(yrs) 


1 


8.3x10'* 


0.10 


1.67 


48.72 




0.015 


0.02 


0.0015 


3 


2.5x10" 


0.01 


0.084 


48.11 


1.0x10"" 


0.91 


0.32 


0.030 


5 


4.1x10" 


0.0037 


0.021 


47.82 


3.1x10*0 


6.03 


1.34 


0.10 


9 


7.5x10" 


0.0010 


0.0043 


47.50 


8.4x1039 


53.8 


6.48 


0.60 


20 


1.6x10'** 


0.0002 


0.0006 


47.05 


1.3x10^' 


894 


50.2 


4.47 



Note. - The HST-1 spectrum used is that of fig.l4l 

"The radius comes from the doubling time (Tablelsl. 
''The standard equipartition magnetic field strength. 

'^R is the ratio of energy densities in synchrotron photons to that of the 
magnetic field. Hence also the ratio of synchrotron self-Compton losses to 
synchrotron losses. 

''The total energy in particles and fields required to explain the observed 
synchrotron emission. 

■^The power required to produce Etot in time dt' (Tablel4l. 

'The halflife in the jet frame for electrons responsible for the UV emission. 
t' is a factor of ten smaller, and these jet frame halflives can be compared to 
the diameter in lightyears to see if geometrical effects or loss times dominate. 

^Tj7d(l.y.) is the ratio of the X-ray half-life to the source diameter measured 
in light years. 

''The halflife in the observer's frame for electrons responsible for the X-rays. 



TABLE 6 
Decay times for HST- 1 



Epoch Interval 


At 


y" 


halftime'' 




(years) 




(years) 


B-C 


0.075 


0.801 


0.19'^ 


G-H 


0.111 


0.817 


0.30 


L-M 


0.126 


0.840 


0.40 


Ya-Yg 


0.115 


0.840 


0.36 


Yh-Yi 


0.126 


0.692 


0.21 



Note. - All values are in the observer's frame, 
"y is the fractional drop of intensity: I2A1 
''halftime=-^ X At 
'■This value was used in Paper I. 



in Table 121 These values may be compared to the estimate of 
Youn2 et al. (2002) for the jet's kinetic energy flow of 3 x 10"*^ 
erg s"' , based on energetics of the interaction between the jet 
and the ISM which are not based on particular jet models. The 
power required to explain the observed flare in M87 would be 
comparable to the jet's kinetic energy only if there were no 
beaming, an untenable hypothesis. 

5. IMPLICATIONS FROM RISE AND FALL TIMESCALES 

There is a long history supporting the notion that bright- 
ness enhancements occur at the locations of internal shocks 
which produce power la w spectra of relativistic ele ctrons via 
Fermi acceleration (e.g.. 'Bla ndford & Konigllll97 9). Polar- 
ization data support this notion with increased polarizations 
and electric vectors perpendicular to the knot leading edge. 
These are characteristic of shocks and ar e associated with the 
flux maximum regions (iPerlman et all 19 99. 2003). 

If we view the shock model as a method of generating a 
new or augmented power law of relativistic electrons (and in- 
creasing the field strength), we note that there is no difficulty 
with the energetics (see above) and also note that the general 
behaviour of the lightcurves in fig. |2] demonstrate compara- 
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ble rates of intensity increase for radio, optical, and X-rays. 
This is a pre-requisite for the shock model which relies on 
the power law extending from low values of 7 up to « 10^. 
Whether the flare owes its genesis to an increase in deliv- 
ered energy from the nucleus, or to a change in the size and 
strength of the shock, is not determined. 

The clear prediction resulting from a new power law distri- 
bution is that the UV emission should decay 10 times more 
slowly than the X-rays, and the radio emission should persist 
much longer because of the 7"' dependence of the electron 
half-life suffering losses. This should be observed in the 
lightcurves so long as the light travel time across the source is 
significantly less than the half-life. For the synchrotron model 
described in Table |5l the estimated size (in l.y.) of the emit- 
ting region is comparable to the halflife of electrons respon- 
sible for the X-ray emission, but less than the predicted half- 
life for UV emitting electrons. In Paper II, we suggested that 
the optical decay times were similar to those in the X-rays, 
but this was based on a quite small intensity drop. Given the 
marked X-ray decline starting in 2005Apr, this issue should 
be clarified by optical and radio data obtained in 2005 but not 
yet available to us. 

6. CONCLUSIONS 

The three lightcurves (X-rays, UV and radio; fig.|3 are par- 
allel to a good accuracy. This means that the input spectrum 
of relativistic particles which caused the flare did not change 
its shape, only its amplitude. The "choppiness" in the X-ray 
light curve may be caused by synchrotron losses which would 
be much weaker at lower energies. 

If we match the observed decay time to the halflife of elec- 
trons producing the X-rays, we find S a bit larger than 5. How- 
ever, if the decay time is governed by the source size, as sug- 
gested by the near equivalence of rise and decay times, we 
should choose a S which satisfies the condition t' < diam- 
eter(l.y.). This requires S < 5, and that would be consistent 
with our results of Paper I that 5 « 4. 

Radio and optical data now being obtained should provide 
the basis for 'solving' two distinct issues: separating light 
travel time from electron halflives and distinguishing between 
compression/expansion versus the injection of new particles. 

If the decay segments of the lightcurves show longer de- 
cay times for lower frequencies, it will be possible to sepa- 
rate the light travel time across the source from loss times. 



From the models already considered, probably the optical, but 
certainly the radio, halflives should be considerably longer 
than the light travel time, so the corresponding decay times 
in these lightcurves will provide a reasonably direct estimate 
of the halflife against E^ losses. Since the magnetic field en- 
ergy density is greater than the photon energy density in all 
models considered, we will thus obtain another estimate of 
the magnetic field strength from the observed halflife. 

If however, all bands rise and fall together, the compres- 
sion/expansion explanation will be supported. Given an emit- 
ting region, be it the result of a shock or some other process, 
a rather modest compression can increase the intensity by the 
required factor of 50 since all electrons will receive the same 
percentage energy boost and the field strength will increase. 
Thus the flare would result from a simple compression with- 
out the introduction of 'new particles'. If the subsequent de- 
cay of intensity were to arise from an expansion, all bands 
should drop equally with no delays between bands. More- 
over, the relative rate of expansion will be reflected in the light 
curve. If the decay timescale is comparable to the rise time, 
this would indicate a rapid expansion (i.e. rapid drop in vol- 
ume emissivity) whereas a decay significantly slower than the 
rise would point to a slower expansion. In the latter case, ex- 
pansion velocities can be estimated from the lightcurve with 
simple assumptions. 



Some of the presented data were based on observations 
made with the NASA/ESA Hubble Space Telescope, obtained 
at the Space Telescope Science Institute, which is operated 
by the Association of Universities for Research in Astron- 
omy, Inc., under NASA contract NAS 5-26555. Work at 
SAO was supported by NASA contract NAS8-39073 and 
grants G02-3144X, G03-4124A, and G04-5131X. Work at 
U. Maryland was supported by NASA through LTSA grant 
NAG5 13065 to the University of Maryland and at UMBC, 
the support came from NASA LTSA contracts NAG5-9997 
and NNG05-GD63G, as well as HST grants GO-9705.01 and 
GO-09474.08. We thank O. Stohlman for writing scripts and 
organizing our flux measurements during the early parts of 
this project. M. Nowak and J. Davis provided useful advice on 
ACIS pileup and the anonymous referee is thanked for helpful 
comments. 
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